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Abstract. The simplicity and the efficiency of a quasi-analytical method for 
solving nonlinear ordinary differential equations (ODE), is illustrated on the study 
of anharmonic oscillators (AO) with a potential V {x) = fix'^ + x^™ (m > 0). 
The method [Nucl. Phys. B801, 296 (2008)], applies a priori to any ODE with 
two-point boundaries (one being located at infinity), the solution of which has 
singularities in the complex plane of the independent variable x. A conformal 
mapping of a suitably chosen angular sector of the complex plane of x upon the 
unit disc centered at the origin makes convergent the transformed Taylor series 
of the generic solution so that the boundary condition at infinity can be easily 
imposed. In principle, this constraint, when applied on the logarithmic-derivative 
of the wave function, determines the eigenvalues to an arbitrary level of accuracy. 
In practice, for /3 > or slightly negative, the accuracy of the results obtained 
is astonishingly large with regards to the modest computing power used. It is 
explained why the efficiency of the method decreases as /3 is more and more 
negative. Various aspects of the method and comparisons with some seemingly 
similar methods, based also on expressing the solution as a Taylor series, are 
shortly reviewed, presented and discussed. 
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1. Introduction 



The accurate determination of the spectrum of the anharmonic oscillator (AO) is an 
old problem which has attracted much interest (for recent reviews see [1]). This 
is certainly due to the asymptotic character of its perturbative expansion [2] (a 
property shared with field theory), but also because it is intrinsically hard to solve 
accurately for sizeable values of the coupling parameter. Of course, the eigenvalues 
may be numerically determined using shooting or relaxation (or similar) methods, 
and extremely accurate results may be obtained that way with the recourse to large 
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computers and sophisticated algorithms (e. g., see [3]). Beyond the fact that a 
purely numerical approach does not allow to learn much on the general properties 
of a solution, it often requires some specific ability and experience. Now in modern 
physics we are more and more confronted with the problem of solving differential 
equations when we are not experts in numerical treatments. It is thus interesting to 
develop simple and efficient analytical (or quasi-analytical) methods which may be 
implemented using symbolic calculation softwares and modest computation powers 

It is not by accident that the mapping method [S] discussed in the present 
paper has been elaborated in the context of the derivative expansion of the exact 
renormalization group equation (for reviews and recent pedagogical introductions see 
[6]). In that field, one is faced with coupled nonlinear differential equations and quasi- 
analytical methods have soon been introduced [71 [8j to solve them because they are 
much more user-friendly than purely numerical methods. Unfortunately they were 
too simplistic and do not converge [21 HHj ■ 

More sophisticated methods have recently been proposed [TTl [5] . They are all 
based on a Taylor series in the independent variable x of the solution of the ODE. 
Methods of that kind had previously been developped in different contexts [lOl [12] . 
Among all these methods, the mapping method seems to be particularly efficient |5]. 
Considering first the logarithmic derivative of the wave function (which softens the 
singularity at infinity), the mapping method relies upon a conformal mapping onto 
the unit disc of an angular sector of the complex plane of the independent variable |5] . 
The angular sector is chosen so that the solution looked for is analytic in its interior 
or assumed to be so. The mapping compactifies the originally infinite integration 
domain whereas the convergent character of the series expansion is transferred to the 
whole disc of unit radius in the new variable. The condition at infinity may thus be 
easily imposed. The object of the present paper is to illustrate the simplicity and the 
efficiency (and also the limitations) of the mapping method with the computation of 
the spectrum of the AO in different configurations. 

The paper is organized as follows. In section[21 the AO is introduced together with 
a change of function (a logarithmic derivative) which transforms the originally linear 
ODE into a nonlinear ODE. In doing so the essential singularity at infinity of the wave 
functions is softened but, at the same time, the new functions are no longer analytic in 
the complex x-plane. The common characteristics of the Taylor-series-based methods 
are then presented. (By the way, the existence of quasi-exact solutions of the AO 
is also shortly illustrated.) Section [3] presents a short review of the Taylor-series- 
based methods utilized by the past to solve the AO. Among those is the popular Hill 
determinant method (adapted to the AO in [13]). Although it is usually not classified 
that way but rather among the variational methods, it is presented here because of 
its close similarity with a simplistic method used sometimes to approximately solve 
nonlinear ODEs (71 [T3]. The two other methods applied to the study of the linear ODE 
are: the power-series method [T^, and the analytic-continuation method [15l[T6]. Then 
methods with an attempt to take explicitly into account the condition at the infinite 
boundary have been used: the Fade method [10] and the contour integral method |12|. 
(The interesting method developed in [17j which is also based on series expansion will 
not be mentioned further because its generalization to any nonlinear ODE does not 
seem possible.) Those latter methods consider the nonlinear form of the ODE which 
induces a non-essential singularity at infinity. In section|4]the principle of the mapping 
method is first briefly reminded. After the transformation of the original linear ODE 
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into a nonlinear ODE, several configurations of tlie AO are explicitly treated with the 
mapping method in view to illustrate its use in different cases. For (3=1 and m = 2 or 
P = and m = 2,3, highly accurate estimates of the eigenvalues are easily obtained. 
For double-well configurations, the accuracy decreases as (3 is taken more and more 
negative (the number of terms in the Taylor series being fixed). The phenomenon is 
related to the decreasing of the radius of convergence of the Taylor series in powers of 
X of the logarithmic-derivative of the wave function. However the method still works, 
in practice more and more terms in the series should be considered to improve the 
estimates when (3 takes on larger and larger negative values. A conclusion is presented 
in section [5] Some numerical results are displayed in the appendix. 

2. The ODE of the anharmonic oscillator 

2.1. The linear ODE 

The usual eigenvalue problem associated with the AO is to find the infinite set of 
values En which correspond to the solutions ipn {x) of the following linear ODE: 

ij" (x) + {E - V (x)) tJj (x) = (1) 

with the condition that the tpn (x) vanish at infinity so as to be squared integrable. In 
([1]) a prime denotes a derivative with respect to x and the potential V (x) considered 
in the present study is: 

V (x) = I3x^ + Ax^™ (2) 

with m = 2, 3, and (3 a given real number. 

With the potential ([2|) , the general solution of H]) is analytic on the whole complex 
x-plane with an irregular singular point located at infinity. 

In the following A can be set equal to unity without loss of generality since the 
energy level of order n, E"™ {/3, A), satisfies the property: 

S,^ (/?, A) = Ai/("+i)£;;" (/3A-2/(™+i), 1^ 

The two-boundary conditions on the corresponding wave functions ipn (x) are: 

la Ipn (x) are either odd or even functions of the independent variable x, i. e. either 
^n (0) = or (0) = 0. 

lb Ipn {x) ^ as a; ^ cx). 

For given f3 and m, the general solution involves three arbitrary constants: two 
integration constants, which are fixed by the parity condition (la) and a global 
free normalization (as it must for an homogeneous equation), and the yet unknown 
"energy" parameter E. The infinite set of values £"„ [the spectrum of llj] is then 
determined by imposing the condition (lb). Actually, without this condition and for 
large x, the general solution of (H]), satisfying the condition (la), has the following 
form: 

^{x) ~ e-^^+Be+^^ (3) 

where B is an arbitrary constant (depending on E) . 

As a matter of fact, determining the spectrum of ([T]) amounts to finding the 
(infinite) discret set of values En (n = 0, 1, • • • , oo) for which B vanishes. 
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Figure 1. Zeroes in the complex x-plane of an approximation of the wave 
function -00 {^) of the fundamental state when /3 = 1, m = 2 with E = 
1.39235164153029 (close to the true value Eq, see the appendix) and for two 
values of the order M of the truncated Taylor series [eq. I I13I I with x = \fz\. 
The approximation of i\> (x) is obtained by simply summing that series. Crosses 
correspond to M = 40 and open circles to M = 80. One sees that the future 
(stable) zeroes of the "true" wave function -00 (a;), are all probably isolated and 
located on the imaginary axis. The other zeroes are potentially not properties 
of the "true" fundamental state, they probably go to infinity as the order M 
increases and as the parameter E approaches the true value Eq. One observes 
among them the possible existence of a positive real zero (which thus goes to oo 
in the same limit). 



Equation ([3]) shows that when S < the general solution vanishes at a positive 
large real value of x. Hence, values oi E ^ £"„ exist for which -0 (x) vanishes on the 
positive part of the real x-axis. 

Figure [1] illustrates the above considerations. 

2.2. The nonlinear ODE 

In order to remove the essential singularity at infinity displayed in ([3]), it is worth 
considering the following change of function: 



in which e = or e = 1 according to the parity of the solution of ([T]) looked for. Then 
the original ODE ^ becomes equivalent to: 



{x) = x'^e^(^) 



(4) 




(5) 



X 
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with the boundary conditions: 



X — >-\-oo 



X 



.m+1 



W{0) = 0, 



W{x) 



m + 1 



(6) 



If instead of W (x) one considers the derivative: 



h (x) = W' (x) 



(7) 



then ([5]) takes the form of a Ricatti ODE (a nonhnear first order ODE). More 
important is the fact that: 



so that any zero of tp (x) in the complex a;-plane becomes a pole for h (x). For example, 
the nodes of the excited states ipn {n > 2) of the AO become poles for h (x) located 
on the positive real a:-axis. 

The zeroes displayed in figure ^ become also poles for h (x) . One sees that 
in the circumstances the true hg (x) would be a meromorphic function of x since it 
would have only isolated poles located on the imaginary axis of the complex x-plane. 
Among the other poles, some are present due to the approximation used (a truncated 
series). But other poles are associated with the existence of a moving singularity 
[characteristic of the general solution of ^] located on the positive real part of the 
a;-axis and which, according to (O [71) is of the form: 



where xo is an arbitrary constant (the value of which depends on the initial conditions 
and of E). 

Determining the true value Eq of the energy of the fundamental state amounts 
to pushing this moving singularity toward infinity. Of course, the location of this 
singularity is directly related to the arbitrary constant B oi 

Via a logarithmic-derivative-like transform ([8]), one is face to the problem of 
solving a relatively simple nonlinear ODE the analytic properties of which are well 
controlled by those of tp (x) solution of the linear second order ODE ^ . This gives the 
opportunity of studying how the various quasi-analytical methods for solving nonlinear 
ODE (presented below) actually work. In general, the estimates of the two first energy 
states Eq and Ei of ([1]) obtained after a logarithmic-derivative-like transform are more 
accurate than that obtained from a direct study of ([T]). This is, presumably, due to 
the fact that the essential singularity located at infinity has been soften. 

2.3. Generalities on the Taylor-series-based methods 

The methods for solving ODE based on expressing the generic solution (i.e., expressed 
in terms of the unknown E) as a Taylor series rely upon the following considerations. 

Consider the ODE ^ as an illustration. It is convenient to perform the change 
of variable x —^ z = x'^ and to redefine the wave function as: 



in which e = or e = 1 according to the parity of the solution of ([T]) looked for. Then 
the original ODE ([!]) becomes equivalent to: 



h (x) = ijj' (x) /ip (x) — e/x 



(8) 




(9) 
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4 z f"{z) + 2 (1 + 2 e) f{z) + [e - V{z)) f{z) - (10) 

with: 

V{z) ^Pz + z"' 

and the boundary conditions: 

/(O) =1 (11) 
/(oo)=0 (12) 

Express / (z) as a truncated Taylor series about the origin z = 0: 

M 

/M(z) = ^a.z^ (13) 



i=0 



in which ao = 1 in agreement with pT|) . 

For foo {z) to be a generic solution of (|10p . the coefficients must satisfy the 
following recurrence relation [obtained by identifying to zero the coefficient of the 
order z'^ in the expansion of equation ^TU\i ]: 

2{k + l)[2k + l + 2e] ak+i + Eak - Pak-i - at-m = (14) 

with the writing convention that = if « < 0. 

Using the recurrence relation p4|l one easily determines iteratively the coefficients 
Qi (E) so that foo (z) is a solution of (fTO|) for arbitrary E. The coefficicents ai [E) so 
determined are polynomials in E. 

Similar considerations stand for the nonlinear ODE ^ and the expansion: 

W{x) ^g{x^) (15) 

M 

gM{z)=Y.c,z' (16) 



i=0 



in which cq = in agreement with ([6]). The corresponding recurrence relation is: 
2fc(2fc - l + 2e)cfc+4^(i + l)(fc-i - l)c,+iCfc_,_i = /34,2--E4,i+4,m+i (17) 

4=0 

in which Sij is the Kronecker delta symbol. From this relation one can determine 
iteratively the coefficients Cfc as polynomial functions of E. 

The expansion of the function h {x) corresponding to the logarithmic-derivative- 
like transform ([5]) is obtained from ([TBI [T7)l using: 

hM {x) = 2xg'M (z) 

with z = a;^. 

To determine the set of values En which potentially correspond to the wave 
functions ijin {x) (satisfying the boundary condition lb) , one has to find one auxiliary 
condition on the coefficients ai (E) or a (E) . The various methods described in sections 
[3] and m differ by the choice of the auxiliary condition. 
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2.4- Quasi- exact configurations 

Let stray from the point a little to stress on the existence of quasi-exact configurations. 
Actually, for certain values of the parameters {/3, m) a part of the spectrum of the AO 
may be exactly determined [18] . Those exact solutions express as polynomials in the 
independent variable x times an exponential. The existence of such solutions has been 
discovered using Taylor-series-based method (Singh et al [19], see also [2Q]). It may 
be interesting to show the mechanism on the particular case m — Z. 

If one factorizes the large x behavior of the decreasing part of ^ before expanding 
in powers of x, then ([9]) is replaced by: 

V'lx) = a;'e-T/(x2) (18) 
and (jl4|) becomes a three-terms recurrence relation: 



2 (1 + fc) (2 A: + 1 + 2 e) hk+i + Sfc - (4 /c + 2 e - 1 + /3) 5fc-i = (19) 

Then, if two successive coefficients vanish, the series terminates to give a 
polynomial. The corresponding solutions are exactly determined. One easily obtains, 
from (fT9|) . the following exact configurations for e = 0: 

/3 = -3, £^ = 0, /(z) = l 

/3 = -7, E^ ±2\/2, /± (z) = 1 T \/2z 

S = ±8, /± (z) = 1 =f4z + 2z2, 



/3 = -ll, 



/(z)-l-fz2 



and for e = 1: 



/? = - 5, S = 0, / (z) = 1 

/? = - 9, E = ±2V&, /± (z) = 1 T 



2^ 



There is one exact configuration for each value of (3 for which the coefficient of 
flfe-i in ([T9)) vanishes (provided fc > 1), namely for the values: 

/3j =-(4j + 3 + 2e), j = 0, 1, 2, 3, ••• 

For each of these values, there are j + 1 exactly soluble eigenvalues whereas the 
remaining eigenvalues of the complete spectrum cannot be exactly determined. In 
section 14.5.31 the quasi-exact configuration (3 = —11 is explicitly treated with the 
mapping method. 



3. The Taylor-series- based methods used to solve the AO 

3.1. The power-series method 

The analyticity properties of the solutions of suggest a very simple method for 
solving the eigenvalue problem of ([T]), namely the power-series method |15j . 
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Using the property that the Taylor series foo (z) has an infinite radius of 
convergence, one may determine the infinite set of values of E for which f^o (zq) 
vanishes at a given finite zq by simply imposing that the sum of the series (jl3p vanishes 
at this point. This zero located on the positive part of the real axis of z is associated 
with a general solution of pO| corresponding to the asymptotic form ([3]) with B 0. 
Then, the limit zq oo should correspond to i? ^ 0. This provides the possibility 
of approaching the spectrum En of the AO (as actually observed in several studies of 
"bounded" oscillators [2T]). 

Dealing with finite values of M , /m (z) represents correctly foe (z) only in a finite 
range of values of z whereas for larger values the sum goes to ±oo according to the 
sign of the last term. The procedure for determining the spectrum of the AO may 
then be sketched as follows. 

2a One fixes M, chooses a consequent but not too large value of zq, and determines 
the zeroes of the polynomial in E associated to the solution of the equation 

/a/ (zo) = (20) 

(this is the effective auxiliary condition at zq). 
2b On increasing zq one observes a convergence of the zeroes towards definite values 

up to a value Zp*^'* where the convergence terminates. At this point one gets the 

best values of the spectrum E„ (z^*''^) for the given M. 
2c One increases M and repeats the process from step 2a. 

2d One observes that Zq*^' is an increasing function of M and that £"„ (^^^o*^'^ quickly 
approaches the spectrum £"„ as M grows. 

This method has been applied in 1961 by Secrest et al T5] with M ^ 1000. They 
notably have determined the fundamental energy Eq with 12 significant figures in the 
case m = 2, /3 = 0. This method has been tried again during the present work, for 
m = 2, l3 — 1 and AI — 250 it gives the fundamental energy Eq with 41 significant 
figures! 

As illustrated by figure ^ , the determination of all the real zeroes of the condition 
(|20p provides an estimation of the complete spectrum, the number of states being only 
limited by the order of the polynomial in E corresponding to the value of M chosen. 
The excited states may be hierarchically determined that way provided zq is chosen 
larger than the location of the last node. As a consequence the number of significant 
figures obtained slowly decreases as the quantum number n increases (for a given Af ). 

It may occur (for some values of the potential parameters (3 and to) that too large 
values of zq or M be required to get a satisfactory accuracy. In that case Secrest et 
al [15j have suggested to reach the value zq in more than one step using the Taylor 
expansion about a non-zero value z/ < zq. This suggestion is already the analytic- 
continuation method introduced later on by Holubec and Stauffer |16] . 

3.2. The analytic- continuation method 

Let us consider the Taylor expansion of / (z) about an arbitrary point z; 7^ 0: 



M 

fM,i {z) = '^h (z - ziY 

i=0 



(21) 
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Figure 2. Part of the even spectrum obtained using tlie power series method 
for /3 = 1, m = 2 and with M = 100. The open circles correspond the real zeroes 
of the polynomial condition ||20|I for varying 20 . The horizontal lines indicate the 
"true" values of the spectrum (Eq, E2, E4, ■ ■ ■ ) as displayed in the appendix. For 
sufficiently high values of zq the points deviate from the lines because fm (2) no 
longer represents correctly the function /oo (z). One sees also that the larger the 
quantum number n is, the larger the minimal value of zq required to identify the 
existence of the energy state must be. 



For foo,i {z) to be a generic solution of (fTOll . the coefficients bi must satisfy the 
following recurrence relation: 

Azi{k + 2) {k + 1) bk+2 = 2 (fc + 1) [2fc + 1 + 2 e] bk+i 

+ {E- Pzi - zD bk-{P + mzY'-^) bk-i 

m-2 , 



E 



bff — m+i 



(22) 



which determines them iteratively in terms of E and of the two arbitrary coefficients 
60 and 61. Those two quantities may be calculated using the series at the origin ([Tc 
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so as to account for the conditions la at the origin with thus: 



^0 = /co {Zl) , 



(23) 



For a given M , having chosen, for example, zi = z — zi — h to he small, one may 
estimate /m (h) and /^^ (h) that way and then /m (2ft.) as a function of the unique 
unknown parameter E using (|21f[23)) . One uses that result to estimate again bo and 
bi at this point in order to reach the point 3h and so on until the point zq — N h is 
reached. 

The following of the procedure is identical to the preceding one (steps 2a~2d 
above) but is relatively slow. Nevertheless it is very efficient provided a well balanced 
choice of M, N and h is done. For example, with M = N = 40, 22 significant figures 
have been obtained this way [22] in the case of the fundamental energy of the double- 
well potential /? = — 50 and m = 2 for which it is almost degenerated with the first 
excited state (the splitting occurs at the 21** figure only). In the circumstances, it is 
more efficient than the mapping method (see section I4.5.3P . 

Notice that, contrary to the power-series method, the analytic-continuation 
method could apply even in configurations where the range of analyticity of ip about 
the origin is limited. 

3.3. The Hill determinant method 

One of the most popular methods used to solve the eigenvalue problem of the AO is the 
so-called Hill determinant method [13] which is (partly) based on a Taylor expansion 
of the solution in powers of the independent variable x such as ((9l [T3|) . 

Actually it is traditionally considered not as a pure power-series-based method but 
rather as a variational method because of an exponential prefactor usually introduced 
in the relation ([9]) between ?/; and the power series (fT3|) . In [T^, Biswas et al use 
a fixed exponential factor e~^ 1"^ but subsequent studies considered e~'^^ with 7 
adjustable [23][24l[25] and even e~''^ with 7 and p adjustable [19]. Nevertheless 
this method may well be sketchily introduced without considering any prefactor. That 
will evidence a similitude with a simplistic method (see section used sometimes 
[7] to solve nonlinear ODEs with the help of power series. 

Let us consider the recurrence relation (|14p as an infinite system of linear algebraic 
equations for the coefficients a^. For this homogeneous linear system, to have a 
solution, its (infinite) determinant (name the Hill determinant, see e.g. [26]), if it 
converges, must vanish. This condition is a transcendental equation for E the infinite 
number of solutions of which should coincide with the complete spectrum of the AO 
provided the condition at infinity is satisfied. 

Of course, practically, one deals with finite values of M and the effective auxiliary 
condition reduces to: 



where Dm (E) is the Hill determinant truncated at order M. 

In fact, the auxiliary condition (|24)) is not always sufficient to determine the 
solution looked for. This is due to the analyticity of the general solution of ([T]): 
the blowing part of the solution also corresponds to a convergent series and may be 
selected by the iterative procedure so defined. In practice one should verify that the 
wave function selected by the Hill criterium (|M]) actually vanishes at infinity. 



Dm (E) ^ 



(24) 
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Using an exponentially decreasing prefactor one may improve the method. It 
is a matter of fact that this "genuine" Hill determinant method works in certain 
circumstances (see for example [IS]). But it fails in certain other circumstances. 
Especially it is unable to furnish the full spectrum of the quasi-exact configurations 
[27] like those presented in section [2T4l see [28] for a review of the reasons of failure. 

3.4- A simplistic method 

At this stage, it is interesting to realize that: 



where du+is) is the {M + 1, 0) minor of the matrix of the system of linear equations 
(|14p . Hence the condition ([M]) is equivalent to the condition that the next coefficient 
aM+i{E) vanishes (e.g. see [191 HI])- This is the condition imposed in certain 
studies of nonlinear ODEs, thus assuming implicitly that the Taylor series converges. 
Unfortunately, in general, this is not true [51 dl] , due to the presence of singularities 
in the complex x-plane [3^ in the case of a nonlinear ODE. 

For example, the solutions of the nonlinear ODE ^ display singularities which 
limit the convergence of the Taylor series ((T6]) . In the case corresponding to the 
logarithmic derivative of the wave function the zeroes of which are displayed in figure 
d] (/3 = 1, TO = 2), one may estimate the radius of convergence of its Taylor series 
from the knowledge of the location of the singularity the closest to the origin (which 
is located on the imaginary axis of the complex x-plane, see figure |T]). In terms of 
z = one obtains a radius of convergece \z\ < i?o — 5.192695. Now by imposing the 
auxiliary condition: 



one assumes that the series converges. But this is only true if [zj < i?o. 

The condition (|25p is not sufhcient to impose the required condition at infinity, 
then it does not yield a convergent procedure for subsequent values of M . Nevertheless 
it, sometimes, gives approximate estimates of the true spectrum. This is because the 
degree of accuracy on the determination of the spectrum is associated with the location 
of the moving singularity. If a large accuracy on Eq is required to push the moving 
singularity beyond the range of convergence of the series then one may expect that (^SJ) 
will give some (limited) information on Eq. In the contrary case, the procedure simply 
does not work. This phenomenon may explain why the simplistic method "works" in 
some case and does not in another case despite similar radius of convergence for the 
respective Taylor series [3T] . 

This simplistic method, when it works, is extremely easy to implement and may 
serve to get a primary estimate of the spectrum before using more sophisticated 
methods such as the following ones which are attempts to construct an auxiliary 
condition that accounts explicitly for the asymptotic behavior of the solution looked 
for. To this end, the singularity located at infinity must be soften, hence these methods 
apply to the nonlinear ODE ([5]). 

3.5. The Fade method 



ttM+l (E) = (-1) 



M+l 




CM (E) = 



(25) 



The Fade method (originally proposed in [10] ) relies upon an attempt to represent 
the solution of ([5]) looked for by successive rational functions of the form: 
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^m,N2 (z) - —pf 

which involves A^i + + 1 coefScients and z — x'^. The processus of the method may 
be described as foUows. 

As is known, the coefficients pi and qi may be determined from the Taylor 
series (fT6| at order M = Ni + N2 according to the usual rules of construction of 
a Pade approximant. This standard procedure determines the coefficients pi and qi 
as functions of the unknown parameter E. The auxiliary condition is then obtained 
by imposing that the Pade approximant so constructed at order M — Ni + N2 still 
reproduces the truncated function at next order M + 1. Namely, setting Ni = s + lo 
and — one gets the following linear system of equations for the coefficients Pi 
and qi (using the convention that c^. — for fc < 0): 

Pi = "^Ci^-jq-j for i = l,---,s + w 
3=0 

s 

= Ci-jqj for i = s + uj + l,---,2s + u! + l 

3=0 

The second line is an homogeneous system of linear algebraic equations for the 
s + 1 coefficients qj . Then the condition to have a solution is that the matrix: 

Tij = Cs+^+i+i_j (i = 0, • • • , s; j = 0, • • • , s) 

must have a vanishing determinant. In the circumstances, Tij is a Toeplitz matrix 
which may be rewritten under the form of a Hankel matrix Hi k by a redefinition of 
the indice k = s — j with thus: 

Hi.k = c^+i+i+k (i = 0, • • • , s; /c = 0, • • • , s) 

Finally, the condition 

det{H) = det ) = (26) 

is the auxiliary condition looked for to determine the spectrum provided that one 
chooses LU in agreement with the boundary condition at infinity. In general it is 
sufficient to choose one of the three values u) = 1,0,-1 according to whether the 
function to be determined goes to ±00, a constant or when z —^ 00. Eventually, 
considering two successive values of uj gives upper and lower bounds on the eigenvalues 

The advantage of the Pade method is that it is well adapted to reproduce, from 
its Taylor expansion, the analytic structure of a meromorphic function especially if it 
has poles. Now, as explain above, the logarithmic-derivative transform changes the 
zeroes of ip {x) into poles for h (x) . Hence the Pade method is better adapted to sum 
the Taylor series of h {x) rather than that of ip {x). Moreover, it is also, a priori, able 
to determine the energies of the excited states despite (or rather because of) the poles 
located on the positive real x-axis (the nodes of i/'n (x) for n > 2). 

The Pade method (named the Ricatti-Pade method in [33] and later on the 
Hankel-Pade method in [35]) has been first introduced in [TU] in conjonction with 
a logarithmic-derivative transform like ^ to calculate, notably, the even and odd 



Conformal mappings for solving ODEs 



13 



fundamental energies of the AO with f3 — 1 and m — 2 and for various values of A. 
Typically the accuracy obtained was about 8 significant figures. Some time later, in 
|32| . the excited state energies have been estimated for the pure quartic and sextic 
AO (i.e. /3 = 0, m = 2 and 3), eleven significant figures were obtained on the estimate 
of £^0 for the quartic AO. This method has then been utilized several times (see e.g. 
[33] for a list of references). In [33] the first two eigenvalues of both the quartic AO 
and the double well down to /3 = — 15 are estimated with 20-18 significant figures. In 
addition, Amore and Fernandez [36] have shown that the Pade method may also be 
applied to solve the two-point boundary value problem associated to several nonlinear 
ODEs. 

In general the Pade method is easy to use [because of the simple writing of the 
auxiliary condition (|26|) ]. It also appears to be robust in several occasions. However 
its effectiveness is limited because the (repeated) calculations of determinants of large 
matrices are extremely time consuming. Sometimes, the Pade approximants introduce 
" spurious" poles or zeroes that can perturbate a clear determination of the spectrum 
En- Though it is not as refined as the following methods, the Pade method may be 
extremely useful (see section l4.5.3p . 



3.6. The contour-integral method 

3.6.1. The fundamental state Leonard and Mansfield [121 137] have proposed the 
recourse to a contour integral in the complex x-plane to perform an analytic 
continuation of the Taylor series of W {x) [satisfying ([5])] towards the large x values 
so that the asymptotic behavior ((6|) can be effectively imposed. The method may be 
described as follows. 

Starting with the Taylor series the coefficients of which satisfy the recurrence 
relation ([rf]). one rewrites this series in terms of large s — 1/x: 

M 

~9M (^) - E i (27) 

According to ((6|), one is interested in finding the values of E for which this series, 
for M ^ oo, has a pole of order m + 1 at the origin s = 0: 

^-(-)^^"- (.. + l).>n+i 

To this end, one considers the following integral over a large circle contour C 
around the origin: 



then, according to the Cauchy formula the contribution of the pole (1^51) to this function 
is — (^jy^^iy (^+i) ■ Assuming that all the other singularities of goo (s) are located on the 
left of the imaginary axis in the complex s-plane then their contributions to F^o (cr) 
will be made negligible as cr ^ oo, so that: 

hm F^ (a) = - ] (29) 

o-^oo (to + 1)! (m + 1) 
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The truncated series ^17} is then used to estimate F^c (u), leading to: 



2i'K& 
M 



1=0 



V{2i + l, 

Owing to the Euler Gamma function in the denominator, this series converges 
(the original series had a finite radius of convergence) and may be summed term by 
term to estimate F^o (c) when a becomes large. Then a procedure similar to that 
described in points 2a to 2d in section [3. II may be applied with zq replaced by cro at 
which point the condition (|29p is tentatively imposed. 

If the region of analyticity of cjoc (s) does not correspond to the assumption that 
all the singularities are located on the left of the imaginary axis of s, the convergence 
may be spoiled by irreductible oscillations. Notice that this is the case for /3 = 1, m = 2 
since, according to figure [H the zeroes of i/'o {x) are located right on the imaginary 
axis. Thus Wq{x) = ln('0o(a^)) has singularities on this axis. To circumvent such 
difficulties, Leonard and Mansfield [T^ [37] propose to modify ([50)) by introducing a 
parameter a/, so that: 



M 

^ r(22a/ + 1) 



22— m— 1 



what corresponds to having performed some rotation of the complex s-plane (on the 
left hand side if a/ < 1). 



3.6.2. The excited states Because the excited states ipn {x) for n>2 have nodes on 
the positive real part of the x-axis, the corresponding Wn {x) have singularities on 
the positive real part of the s-axis which cannot be moved by a/. Hence the method 
does not apply directly to the determination of the excited states. In order to have 
access to them, Leonard and Mansfield [T? propose to use Wq [x) (determined by 
the procedure described just above) as a basis to write [here accounting for the odd 
(e = 1) and even (e — 0) possibilities]: 



V'(a;) =x^e'^'(^)p(a;) 



(31) 



The ODE satisfied by P (x) is then: 



P" + —P' + 2W'^P' + {E- E,)P^O 



This differential equation has two types of large x solution: 

E- E, 



Fix) 



P{x) 



exp 



exp 



2(m- l)2:'"~i 



(32) 



(33) 



(m+ 1)_ 

Only (|32p is compatible with the boundary condition at infinity for ^„ {x). 
Actually, because in ([3T|l We [x) already factorizes the exponential decreasing of 
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V'n (x), Pn (x) must approach a constant, then p2p satisfies this condition since m > 1. 
Taking into account one sees that the second behavior ([55)1 reconstructs the 
blowing part proportional to B in ^ as it must. 

With the behavior ([5^ . P„ (1/^) ^las no singularity at s = and the contour 
integral procedure may again be applied to determine the values of En — E^ that make 
this integral vanish when cr ^ oo. 

Notice that the asymptotic condition (zero for the excited states) being 
independent of the energy state, an individual adjustment may correspond to any 
value En so that a global (even partial) obtention of the spectrum [as that illustrated 
by figure ([2])] is not possible. In particular, a plot of the wave function is necessary 
to identify the excited state so determined [12]. One could, however, use a procedure 
similar to that of section 13.11 for the power-series method to determine globally the 
spectrum before looking for higher accuracy on the determination of each individual 
energy state. 

3.6.3. Resettling In order to improve the efficiency of their method applied to the 
AO, Leonard and Mansfield [12l[33 utilize a rescaling which allows them to assign to 
the third coefficent C3 the role of the adjustement parameter instead of E whereas C2 
is fixed instead of /?. With this trick and 300 terms in the series (Af — 300), they 
are able (for C2 fixed to zero, hence (3 = Q and to = 2) to estimate Eq with an 
accuracy of 65 significant figures. The remaining of the spectrum of the AO up to 
n ~ 40 is determined with an accuracy of 48 significant figures. That is a remarkable 
achievement. 

It is to be noted that the rescaling trick is not very convenient when /? 7^ since 
then a supplementary adjustment of C2 is required [12]. 

As indicated below, the mapping method appears to be more efficient than the 
contour-integral method, since, using of the same rescaling trick, it yields a much 
greater accuracy whereas M = 250 only for the configuration /3 = and m — 2 (see 
the appendix). 

The major drawback of the contour integral method is the iterative uncomfortable 
adjustment procedure necessary to approach the infinite boundary (similar to that of 
the power-series method described in section [XT]) . 

4. The mapping method 

4.1. Introduction 

In solving the AO problem, it is important to take correctly the condition at infinity 
(an essential singularity) into account. It is thus useful to first reduce this singularity. 
To this end a logarithmic-derivative of the wave function is usually performed. As 
illustrated with the Pade and the contour-integral methods, the condition at infinity 
may then be more efficiently imposed. With the contour integral method an analytic 
continuation of the original Taylor series towards the infinite boundary is used. But 
no attempt is made to compactify the integration domain of the original ODE, so that 
the practical execution is finally similar to that of the analytic continuation method 
(the power-series method even) with a clumsy step-by-step approach to a point located 
at infinity. 

The mapping method introduced in [5 leans also on an analytic continuation after 
a logarithmic-derivative-like transformation but in addition, the infinite boundary is 




Figure 3. When the interior of an angular sector involving the positive real axis 
of z defines a region of analyticity for the function / (2) (with 2 = x^), it may 
be mapped onto the interior of a circle of unit radius for the complex variable w 
defined by eq 1134 l l. The singularity 20 of / (2) the closest to the origin determines 
the analyticity domain provided the other singularities lay outside the angular 
sector containing the origin. 



brought close to the origin using the foUowing conformal transformation: 



1 



(34) 



{1 + z/Rf" 

in which R and a characterize the position of the vertex and the angle of an angular 
sector of the complex plane of z as shown in figure [3l 

The conformal transformation (|34p maps the interior of the angular sector of the 
z-plane into the interior of the unit circle centered at the origin of the w-plane so that 
z = 00 corresponds to w — I (whereas z = corresponds to w = 0). 

If the interior of the angular sector is a region of analyticty of the original 
function g (z) then the series ()16|) converges there. Under the analytic continuation 
this convergence is conveyed to the whole unit disc for the series: 



M 



(35) 



obtained by expanding, within the original series ()16|) . the relation inverse of ([34]): 

'1 + w 



R 



1 — w 



1 



(36) 
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M 

Figure 4. Evolution witli tlie order M of all the real zeroes corresponding to the 
auxiliary condition II37I I for the AO with /3 = 1, m = 2. The mapping transform 
1134 l l is used with R = 5.192695 and a = 2 corresponding to the expected analytic 
structure of the fundametal state i/iq (see figure [T] and text). The candidate to 
the fundamental energy Eg is well evidenced. 

4-2. Example of the AO with (3—1 and m — 2 

Figure [T] shows that, presumably, aU the zeroes of the fundamental state ipo (x) are 
located on the pure imaginary axis of the x-plane. Hence, in the z-plane of the function 
g (z) , which is related to ip {x) via (^1 [T51 [TB|) . those zeroes become singularities located 
on the negative part of the real z-axis. The distance to the origin of the closest 
singularity provides the value Rq ~ 5.192695 and the plane cut on the negative real 
axis, starting from the point z ~ —5.192695, forms an angular sector such that R — Rq 
and a = ao = 2. Choosing those values in the conformal mapping (j341 136p . produces 
a truncated series gM (w) in powers of w. As M is increased the sum of g]\j (w) at 
w = I should approach the asymptotic behavior corresponding to ([6l I15[ fT6|) . Thus, 
accounting for those constraints, the auxiliary condition looked for to estimate Eq may 
be expressed as: 

-g'ii' H^=i = (37) 

in which g^"^^ (w) stands for the conformal mapping applied on the function g" (z) 
which, according to ([51 [HI [H]) , goes to zero as z ^ cx). 

The condition p7p is a polynomial equation for E, the solutions of which 
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M 

Figure 5. Variation with the order M of the approximate number of stabiUzed 
figures (A^c = — log\l — Em /Em— i\ where Em—i and Em are two successive 
estimates of Eq) in the determination of the fundamental energy Eg of the AO 
with /3 = 1, m = 2. The method used is the conformal mapping of (|34^ with 
R = 5.192695 and a = 2. 

effectively display a stable real value as M is increased. Figure ([4]) shows that, even 
when M is small, this value is easily identified and can be followed without difficulty. 
Actually this stable value converges quickly to the true value Eq. As example, with 
M = 150, Eq is determined that way with an accuracy of 83 significant figures (see 
figure [5] and the appendix). Only the choice of the order AI limits the accuracy of 
this estimation which is much more accurate than the previously published values 

[islElllsH]. 

An efficient and very convenient variant to the auxiliary condition (I37p consists 
in simply imposing the vanishing of the last term of the series ([35]), say 

UM iE)^0. (38) 

This condition avoids the explicit account of the asymptotic behavior of ^pQ (x) . 
That is quite justified knowing that the general solution (for a generic E) has a 
moving singularity located right on the positive real w-axis which generally limits 
the convergence of the series gM (w) except for the value Eq for which this moving 
singularity is sent to infinity. Hence imposing the condition ([38]) amounts to force the 
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Figure 6. Illustration, on the AO with /3 = 1, m = 2, of the respective 
importance of the free parameters R and a, the correct values of which would 
be Ro ~ 5.192695 and a = 2. Nc is defined in the caption of figure[5] On the loft, 
a is fixed to 1 and R is varied. It is clear that the best convergence corresponds 
to -R = 5. For smaller values of R, the convergence seems to be monotonically 
increasing whereas for i? = 7, it seems to be limited. On the right, R is fixed to 
5 and a is varied. The best convergence seems to correspond to a = 5/2 which is 
a too large value. Actually, the small rule box shows that at larger values of the 
order M, this choice, most probably, is no longer the best one. 



convergence of the series and this procedure enable the determination of the value of 
Eq. 

4-3. Practical use of the mapping method 

In general one does not know a priori the values Rq and ao which correspond to 
the location of the singularity the closest to the origin of the solution Wq (x) looked 
for. One must thus consider R and a as free parameters. Several procedures may be 
conceived to approximately determine the "best" values of R and a (which, ideally, 
should be close to Rq and ao respectively). The following one is illustrated on the 
case (3=1 and m = 2. 

Suppose first that Rq and ao define an angular sector in which the solution looked 
for is analytic (i.e., there is no supplementary singularity in the interior of this sector). 
Then, if R < Rq and a < ao, the mapping method should provide a convergent result 
as the order M is increased. If it does not, then one of the previous inegalities is 
not satisfied and a decrease of the trial value of R (also of a even) is necessary. One 
thus proceeds by successive trial and error in order to determine at least one couple of 
value {i?i,ai} such that the mapping procedure begins to converge. In the example 
presently chosen, if one sets a priori R = 1 and a = 1, one already gets, with AI = 50 
an estimate of Eq with 8 significant figures (1.39235163). Having this first estimate, 
one may look at the effective radius of convergence of the original series. This may be 
done with the help of the d'Alembert or the Cauchy rule or (in the case where they 
do not converge) via the plot of the sum of the series as function of z. This provides 
an estimate of Rq which may be used to get a better convergence towards Eq. It is 
important to realize that, to get a relatively large accuracy on Eq, it is not necessary 
to know Rq with much accuracy. 

The determination of the best value of a proceeds also by trial and error and is 
determined by the criterium of best convergence as shown by figure However, 
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Figure 7. Evolution with M of the distribution of the real zeros (open circles) 
corresponding to the auxiliary condition of the mapping method applied on the 
function P {x) defined by 13111 in the case /3 = 1 and m = 2 (even case). The 
horizontal lines correspond to the "exact" values of the spectrum given in the 
appendix. 



for small values of M, it may appears that "exotic" effective best values of a may be 
observed. Figure ^ shows that in the example of the AO presently treated, the value 
a — 5/2 provides, for small M , an apparent better convergence than the right value 
a = 2. But, for larger values of A/, one observes that the former case finally yields 
perturbed convergence whereas with the second case the convergence remains smooth 
[see figure ©]■ 



4- 4- Excited states 

The determination of the excited states proceeds as in section [3?6l The fundamental 
state Wq {x) (or the first odd state Wi (x)) is used as a basis through (pij) . 
Consequently, there is no need to look for new determinations of R and a. The 
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values obtained previously for Wo,i (x) work also for the excited states. At a given 
order M , the corresponding auxiliary condition yields a polynomial equation for E, 
the real zeroes of which form, once and for all, the approximate spectrum of the AO 
which can be, step by step, unambiguously determined as shown by figure ([7]). The 
values of the spectrum are given in the appendix, they are much more accurate than 
(and compatible with) the existing previous estimates [13l [211 [38] . 

4-5. Other examples treated 

In this part configurations which are more difficult to treat than the previous one are 
considered. One easily understands that the larger m, the harder the determination 
of the spectrum. In principle the difficulty also grows as (3 decreases: for a fixed m, 
the limit /3 ^ is the limit of infinite coupling constant A ^ oo. Another kind of 
complication arises when f3 is more and more negative (double- well configurations): a 
Taylor-series-based method is not very appropriate because the expansion is performed 
on the top of the hill of the potential. 

4-5.1. The AO with (3 = and m = 2 This is the example treated explicitly in 
[T^ where, with AI = 300, the fundamental energy Eq had been determined with an 
accuracy of 65 significant figures (48 figures for the excited-state energies En up to 
n = 39). 

The analytic structure of Wq (x) is similar to that of the case /3 = 1 previously 
considered: all the singularities are probably located on the imaginary axis of the 
complex a;-plane. The location of the singularity the closest to the origin implies 
i?o ^ 4.2526611666886 whereas ao = 2. 

Using those values to fix R and a, the mapping method yields, for M = 150, an 
estimate of Eq with already 76 significant figures. Moreover, using the rescaling of 
fi2\ [37] and a series calculated up to the order M = 250, the number of significant 
figures obtained is 114. The accuracy is even better in the odd case since, at the same 
order M = 250, Ei is estimated with an accuracy of 120 significant figures (see the 
appendix). 

As for the excited states both in the even and odd cases and M = 180, the 40 
first energy states are estimated with an accuracy larger than 56 figures in the even 
case (63 in the odd case). 

All these estimates are presented in the appendix. They all coincide with the 
estimates given in [I^ but, again, involve more significant figures. 

4-5.2. The AO with P = and m = 3 In that example, the singularity of Wq (x) are 
not located on the imaginary axis of the x-plane but on four straight lines starting 
from the origin so that the singularities the closest to the origin are located at points 
re^''^ and re=^*(''-^' with: 

r ~ 1.81997592362 

6 ~ 0.3663025 x tt 

Considering now the complex z-plane, and the generic definition of R and a 
associated with an angular sector as drawn on figure ([3|), this induces the angular 
sector defined by: 

Ro ~ 3.31231236256 

ao ~ 0.732605 
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with two singularities located on its edges. With such a sector of analyticity, without 
singularity taking place on the negative side of the real 2;-axis, the parameter R is 
no longer bound by the condition R < Rq and may be chosen as large as one wants 
provided that a be modified adequately. However the larger R, the smaller a, and it 
appears that it is preferable to choose the largest integer value of a which is compatible 
with the constraint R > 0. Presently the value a = 1 may be chosen whith R ~ 2.2107. 
With those values and the rescaling of [12l[37] , using M = 250, the fundamental energy 
Eq has been determined with 56 significant figures. The 20 first even excited states are 
then determined with at least 24 significant figures with M = 181. Similar accuracies 
are obtained in the odd case (see the appendix). 

4-5.3. Double well The more f3 is negative, the more the calculations based on a 
Taylor expansion about the origin a; = are difficult. In addition, one will see that the 
radius of convergence of the Taylor series for Wo,i (x) decreases so as to considerably 
reduce the efficiency of the mapping method (although it continues to work correctly) . 
For small (3 < and m = 2, the mapping method is still very efficient as shown by 
the estimates given in the appendix for (3 — —1 and m = 2. 

For the sake of an illustration of the difficulties encountered for increasingly 
negative values of /3, the following is limited to a short presentation of the results 
for two configurations: /3 — —11, m — 3 and (3 = —10, m = 2 (the configuration 
/? = — 50, m = 2 is also shortly alluded to). 

/3 — —11, m — 3 According to section that configuration is quasi-exact. One 
finds that the singularities of Wq (x) are all located on the imaginary axis in the 
complex x-plane and the optimal values of the mapping parameters are determined 
to be: 

Ro ~ 0.29289 

ao = 2 

According to ([6]), the auxiliary condition is: 

g'i^ H^=i = -I (39) 

With M = 100, the estimate obtained with the mapping method for the fundamental- 
state energy is: 

Eo = -7.99999999999999969 ± 1.4 x 10"^^ 

which corresponds to the first exact value —8 listed in section [^T^ 

In the odd case and M = 100 also, Ei (which is not an "exact" value) is better 
determined with an accuracy of 30 significant figures, using Rq ~ 0.9129007 and 
QfQ = 2. It comes: 

El = -7.91735504373511193638182033627 

To get the two other exact values, one must consider the even excited states. 
With M = 100, the results are: 

E2 = 0.00004 ±0.0007 
E4 = 7.99976 ± 0.00397 



Conformal mappings for solving ODEs 



23 



which correspond to the exact values and 8 of section 12.41 respectively. 
The two following odd excited states are again better determined: 

E3 = 2.520358324953 
E5 = 14.11296433804 

To improve the accuracy, one should have to consider larger values of M. 

One notices that, compared to the even case, the better accuracy obtained in 
the odd case coincides with a larger radius Rq. The following example confirms this 
(expected) correlation. 

P — —10, m = 2 The singularities of Wq (x) are again all located on the imaginary 
axis in the complex x-plane and the optimal values of the mapping parameters are 
found to be: 

Ro ~ 0.118685154 

ao =2 

According to ([5]), the auxiliary condition is: 

ffi'^ = (40) 

The radius Ro is smaller than in the previous case and a less accurate estimate 
must be expected. Effectively with M = 126, the estimate obtained with the mapping 
method for the fundamental-state energy displays only 9 significant figures: 

Eo = -20.63357672 

and for M = 200, one gets (after a relatively long time of calculation) about 16 
significant figures only: 

Eo = -20.633576702947807 
In the odd case, one finds: 

Ro ~ 0.44986354 

ao =2 

and the estimate with M = 100 is found with an accuracy of 16 significant figures: 

El = -20.63354688440494 

The above estimates of Eq and Ei are not as accurate as those obtained using 
the Fade method with M ~ 40 only [33]: 

Eo= - 20.633576702947799150 
Ei= - 20.633546884404911079 

One sees that, in accordance with its nature, the Fade method makes the best 
with the (isolated) singularities located on the negative real z-axis even if one of them 
comes close to the origin. On the contrary, the efficiency of the mapping method 
decreases in that case because it is intimately tied to the location of the singularity 
the closest to the origin. 
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P = —50, TO = 2 It is worth mentionning that, in that case, the radius Rq is found 
to be close to 0.004012. With a so small value, the mapping method is practically 
useless (although it continues to be valid in principle but large values of M are already 
required to simply get an idea on the order of magnitude of i?o). In the circumstances 
the analytic-continuation method of section l3.2l is particularly useful since, as already 
mentionned, that method yields an estimate of Eq with 22 significant figures [22^ with 
M ~ 40 only. 

5. Summary and conclusion 

The AO with different potential configurations have been explicitly considered to 
illustrate the efficiency and the limitations of a quasi-analytic method for solving 
nonlinear ODE named the mapping method [5]. As several other quasi-analytic 
methods encountered here and there in the literature of the AO it is based on a generic 
Taylor series in powers of the independent variable. After a short presentation of those 
methods, the mapping method has been introduced and its use clearly illustrated on 
the basis of the simplicity of the analytic properties of the AO. In particular the 
practical determination of the two adjustable parameters R and a, inherent to the 
method [see eqs. (|34l l36| ]. has been shown to correspond precisely to those analytic 
properties. In absence of any information on them, the criterion of best convergence 
is appeared to be valid. Extremely high accurate estimates of the spectrum of the AO 
have been easily obtained (see the appendix) provided that the effective value of R 
was not "too" small since, in accordance with its nature, the efficiency of the method 
decreseases with R. In the case of very small R, the analytic-continuation method or 
the Pade method have appeared to be more efficient. Finally the mapping method 
[S] is an extremely refined method, easy to use and which provides clear convergences 
towards the values looked for provided the radius of convergence of the initial Taylor 
series is not too small. The Pade method [10] is also an easy-to-use method which is 
more robust but less refined than the mapping method. The two methods (and also 
the analytic-continuation method J15','16]) may certainly be advantageously associated 
in the process of solving a two-point boundary problem of a nonlinear ODE. 
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Appendix 

In this appendix are presented the estimates of the energy spectrum of several 
configurations of the AO as obtained using the mapping method. 

Appendix A.l. Spectrum for /3 = 1, to = 2 

With M = 150, R = 5.1926948467766235669 and a = 2, the fundamental energy state 
has been determined with 83 significant figures: 

Eo = 1.39235164153029185565750787660993418460006 
67112208340889063493238775674318756465286 
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n 

2 8.655049957759309688116539457377308026275 40 

4 18.0575574363032528947712396465254348531 39 

6 28.83533845950424884013363571549983817 37 

8 40.690386082106444725278931481582464 35 

10 53.4491021396652646008315064597595 33 

12 66.995030001247166061019704904702 32 

14 81.24350505076715273706652147034 31 

16 96.129642045234052046811222396 29 

18 111.6018150451729585337015116 28 

20 127.617777795354918333962292 27 

22 144.142195296398163731983 24 

24 161.14499069451295186862 23 

26 178.6001923668757611938 22 

28 196.48510291022044366 20 

30 214.779683549176627 18 

32 233.4660874793752 16 

34 252.5282990614935 16 



Table Al. Estimates of even excited state energies fo /3 = 1, m = 2 as obtained 
by the mapping method with M = 115. The last column displays an approximate 
value of the number of significant figures obtained in each case. 

The estimates of the 17 first even excited states, obtained with M — 115, are 
given m table EH 

With M = 115, R = 6.03344983950017 and a = 2 the first odd state has been 
determined with 70 significant figures: 

El = 4.64881270421207753637703291726058448 
8898860447882825934823424910341006 

The estimates of the 20 first odd excited states, obtained with M = 121, are given 
in table El 

All the estimates obtained are in agreement with the existing literature [551 [351 US] 
at least up to 10-16 significant figures they quote. 

Appendix A. 2. Spectrum for /3 — 0, to = 2 

In order to compare the efficiency of the respective methods, the rescaling of [19] 
has been used in the study of this configuration. In the even case, the parameters 
of the conformal mapping have been fixed to i? = 0.7515601141384727 and a — 2 
corresponding to the location of the singularities of Wq {x) in the complex x-plane. 
With M = 250, the fundamental state energy has been determined with 114 significant 
figures: 

Eq = 1.06036209048418289964704601669266354551 
52087285289779332162452416959435630443 
4442112689629913467170351054624435858 

The 20 first even excited state, with M = 180, have been determined with an 
accuracy as displayed in table IA3I 
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n Er, 

3 13.1568038980498750792097720403823146746501484 45 

5 23.2974414512231890848644819920981238281208 42 

7 34.640848321111332542884527618156342033769 41 

9 46.9650095056755279840964433241751142524 39 

11 60.12952295915777131584801605985282216 37 

13 74.03587435910253018074120548740369915 37 

15 88.61034880079915887303910537132488 34 

17 103.7953003222726096781116879557136 34 

19 119.54417073305031113002694934564 32 

21 135.818417325610373340451430114 30 

23 152.5855042055739215668661903 28 

25 169.817528001595348199877321 27 

27 187.4902426929503225448058 25 

29 205.58234660442351871834 23 

31 224.0749478526003062853 22 

33 242.95115495114712353 20 

35 262.1957574685198472 19 

37 281.794972923819312 18 

39 322.008069744848 15 

41 363.5018948643 13 



Table A2. Estimates of odd excited state energies fo /3 = 1, m = 2 as obtained 
by the mapping method with M = 121. The last column displays an approximate 
value of the number of significant figures obtained in each case. 



In the odd case, the mapping parameters are fixed to R — 0.65325380587388 and 
a = 2. The first odd energy state has been determined with 120 significant figures 
(M = 250): 

^1 = 3.7996730298013941687830941885125689577660 
6546732735340998309740435228060023769937 
931149834764800025936072118065726127409 
With M — 181, the 20 first odd excited states are estimated as shown in table 

Ell 

All the above values are in agreement with the estimates published by Leonard 
and Mansfiled [T^] at least up to 48-65 significant figures they quote. 

Appendix A. 3. Spectrum for /3 = 0, to = 3 

In this case the singularities of VFq {x) are not located on the imaginary axis of the 
complex x-plane. Hence a 7^ 2, and R may be adjusted so that one may choose a = 1 
to make the calculations easier. Using the rescaling trick of [19], for i? = 0.42181 and 
M — 250 the fundamental energy has been determined with 56 significant figures: 

Eo = 1.1448024537970527637654575341495490765378125289798777292 

With M — 181 the 20 first excited even energy states are determined to be as in 
table Ell 



Conformal mappings for solving ODEs 



27 



n 




En 


Nc 




7.455697937986738392156591347185767488137819536 


86 


2 




74910086643688870645553641487360130430909 




16.26182601885022593789495443038461353424458650 


86 


4 




45141645421211209355815531735099805375715 




26.52847118368251819181382818368122792897571351 


84 


6 




294412586757107839537851365173102739202 




37.92300102703398514651637855190995445547624176 


82 


8 




6768540940158056680472093643981051811 


10 


50.25625451668291903974458810526275993610124418 


81 




075701487761224602795878984785502928 


12 


63.40304698671889497209348492294463561729801089 


80 




28673988084514263379591656250007033 


14 


77.27320048198396570945437444123301297760205161 


78 




696448894060040064817405480052452 


16 


91.79806680899113406984882666252278502713173398 


77 




90657479115965604376782712914156 


18 


106, 


.9233073817325256530751021361667347590846183 


76 




2008357932027830540845145707859 


20 


122, 


.6046390009994550207629714176151818749766332 


74 




22893621015463272615565590626 


22 


138, 


.8051479113947530940883610374380934368840381 


72 




671106897341615306983140020 


24 


155, 


.4935022686820551465506277045177934315466954 


70 




5414662134066104062468041 


26 


172, 


.6427119628444830668784007249101646087468190 


68 




46656519547351087844600 


28 


190, 


.2292386524627088279661423652564958817125565 


66 




189811610764662196522 


30 


208, 


.2323390051439485027234771412174093538259896 


65 




12647599436328976471 


32 


226, 


.6335684811375306036977950956875757647131141 


64 




4054544586360420784 


34 


245, 


.4163987919363679916747517301130143224691180 


61 




4034148423591243 


36 


264, 


.5659178144988804976115801604901826931113435 


60 




585310082273247 


38 


284, 


.0685905814007431504962812081250647770847132 


58 




6709967957974 


40 


303, 


.9120663483839986502856546215441904999429332 


56 




70794969657 



Table A3. Estimates of even excited state energies fo /3 = 0, m = 2 as obtained 
by the mapping method with M = 180. For each n a second Une displays the 
continuation of digits. The last column displays an approximate value of the 
number of significant figures obtained in each case. 
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n 


En 


Nc 




11.64474551137816202085037328137093643655087216197 


y4 


3 


4453001669270913851625610586869623215876947849 




21 .23837291823594002414971 1 11358863637670483205965 


92 


5 


52359735931435198551551347510844880262694475 




32.09859771096832663427210643833189668755992423305 


90 


7 


978762342171173242882982430096400348694857 




43.98115809728973078531811375282736479477947632092 




9 


7591678745311041539444158927512411142579 


88 


11 


56.73421405517303604730167103345175751218649581511 




39323399711320613883016063024691753961 


86 


13 


70.25239462861659088913085193801913436712112961025 


85 


6667440126160413192685489535915017187 


15 


84.45746627494195892888833727655962244102394532902 


83 


36232438363406557704202821636803725 


17 


99.28860666049327964581980672652356738129618429531 


82 


5439623923490779365455043123425144 


19 


114.6969173849851755381509283197832144905718574275 


81 


094513375030013226382190538593688 


21 


130.6420687486299788272084813607924461800377164913 


80 


39742213272421838252016675620413 


23 


147.0901212576038250593001573616393276076304105930 


78 


039522847501664310955551882108 


25 


164.0120436228653346009732131640724071039525335388 


76 


6631868425810048601661780456 


27 


181.3826661857682257850039425558429144469164137642 


74 


86635610868225155373630509 


29 


199.1799188337470723626128506406160415222128605618 


72 


585703827731174998592784 


31 


217.3842616741034519777291998578486546497746509056 


71 


71607931000643490074315 


33 


235.9782503616957543907796533463549468434869938281 


69 


882674294457757985404 


35 


254.9461979707982589533328548019610045320736538738 


68 


26605804407445584226 


37 


274.2739076589405355230240455122439840737715234577 


66 


772581092942928166 


39 


293.9484582660060854336699974835216263034458992746 


64 


8022674262238544 


41 


313.9580301839779322009488360107644470523106026887 


63 


398971572795076 



Table A4. Estimates of odd excited state energies fo /3 = 0, m = 2 as obtained 
by the mapping method with M = 181. For each n a second hne displays the 
continuation of digits. The last column displays an approximate value of the 
number of significant figures obtained in each case. 
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n iV, 

2 9.073084560921433856016249096664712179552456747 46 

4 21.7141654221967222816897851486065849098555916 45 

6 37.6130865608951608896015430474867021943435 42 

8 56.199300852499359421090060527088020427952 41 

10 77.12734146382543798734243263215304574937 40 

12 100.15892789265287617620112949318711203 38 

14 125.11668074758135609314911179285667525 38 

16 151.86183175856232716716355175727249217 38 

18 180.28193894068557040564022178012003 35 

20 210.2834723116004980559950901608572 34 

22 241.787034056498822149736098165624 33 

24 274.72411811463802668627165975505 32 

26 309.03482436375835367061889961183 32 

28 344.666193379464295659578007932 30 

30 381.57096047135932038713367165 29 

32 419.7066022261534757147058497 28 

34 459.034592731153881282354894 27 

36 499.5198136445789902188947 25 

38 541.130079449256802989435 24 

40 583.835750482220100204957 24 



Table A5. Estimates of even excited state energies fo /3 = 0, m = 3 as obtained 
by the mapping method with M = 181. For each n a second line displays the 
continuation of digits. The last column displays an approximate value of the 
number of significant figures obtained in each case. 

A similar singularity distribution occurs for the odd state so that a may be fixed 
to 1 with R = 0.3467407074263. With M = 250 and the rescaling trick of [19|, the 
first odd energy state is determined with 59 significant figures: 

El = 4.3385987115139811916473368880601867793285474529432549486123 

With M — 182 the 20 first excited odd energy states are determined to be as in 
table Eel 

All the above values are in agreement with the estimates published in [38] at least 
up to 16 significant figures quoted. 

Appendix A. 4- Spectrum for f3 = —1, m = 2 

In this case the zeroes of the two first states (even and odd) are located on the 
imaginary axis of the complex x-plane and the parameter a of the mapping transform 
(|34p is set to 2. The values of R corresponding to the location of the zeroes the closest 
to the origin in the two cases are found to be: 

R = 3.33794474927577803 (even) 

R = 4.31581775651355421 (odd) 

With M — 100, the tables [ATI and lA8l give the values of the 21 first energy states 
in each cases (even and odd): 

These results are in agreement with the estimates given in the literature 
[39l l40l [22l [TBI [38] at least up to the significant figures they quote. 
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n iV, 

3 14.93516963491073601571195495127066852852812677095 49 

5 29.29964593740189312978407064699053036251622441 46 

7 46.59521144855173401972883161023883550190357379 46 

9 66.387281706591605561901323787042869206296808 44 

11 88.392375769030215606669098836449373556533 41 

13 112.40643581144162471251824362883088085074 41 

15 138.2733726559526082048335446213274433132 40 

17 165.86873924695253971482736881612442986 38 

19 195.0902772336198933064947454865236007 37 

21 225.852006093743593994689769811061848 36 

23 258.08030930413349189359807845604983 35 

25 291.7112248219607946631497420954846 34 

27 326.688501175539696874664982050866 33 

29 362.96216142162740255082506466132 32 

31 400.4874160213062718461371957159 31 

33 439.223822596901596341934330594 30 

35 479.13462481132026473711629662 29 

37 520.186224052830671932702626 27 

39 562.347751461637717237705920 27 

41 605.59071704038135107740483 26 



Table A6. Estimates of odd excited state energies fo /3 = 0, m = 3 as obtained 
by the mapping method with M = 182. For each n a second line displays the 
continuation of digits. The last column displays an approximate value of the 
number of significant figures obtained in each case. 
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